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Abstract 


The large signal behavior of linear beam travelling wave 
amplifier is presented in this paper. The electron transit time 
and the electron exit velocity in a dielectric supported helix 
traveling wave tube amplifiers (TWTA) are derived by solving 
the non-homogeneous boundary value problem. From the 
Maxwell’s field equations, electron ballistic equation and from 
the boundary conditions of the tape helix TWT, a Green’s 
function sequence is obtained for the axial electric field as a 
function of electron arrival time in the interaction region. By 
Successive approximation technique, the electron transit time 
and therefore the exit velocity equations were derived and 
numerically computed. 


Keywords: Electron ballistic equation;  anisotropically 
conducting, tape-helix, electron beam; boundary conditions 
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Introduction 
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Fig. 1. Cross sectional view of Dielectric supported Tape Helix Travelling Wave Tube (TWT) 


Applications of TWIT: High power satellite transponder output, 
medium power satellites, final amplifier for radar system. 


All structures are redrawn based on the structures present in www.cst.com 
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Working Principle of Helical TWT 
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Fig. 2. Cross sectional (y-z plane) view of interaction region of the slow wave structure 
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Electron Beam Interaction with r.f. fields 


> DC beam velocity of the beam is maintained slightly greater than that 
of the axial field. 

> The helical slow wave structure (SWS) helps in the interaction of the 
electric fields and the electron beam in the interaction region, by slowing 
down the velocity of the RF fields. 
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Fig. 3. Electron bunching due to accelerating and decelerating fields 1n interaction region 
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Analysis models of tape helix conductors 
Sheath helix model: Tape width, w = p, pitch length. Current direction is 
assumed to be along the tape direction. All existing large signal analysis of 


slow wave structures involve sheath helix model. 


Tape helix model: Tape width, w < p, pitch length. This is more practical 
approach. 
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Fig. 4a. Sheath Helix Model Fig 4b. Tape Helix Model 
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Tape helix model: Anisotropically 
conducting tape helix 





Anisotropically conducting tape helix model has surface current density only 
along the direction parallel to the tape winding direction. (Jy) 
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Fig. 5. Anisotropically conducting tape helix structure 
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Azimuthally averaging the region 
between tape helix and outer conductor 
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Fig. 6. Cross sectional (x-y plane) view of interaction length 
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Fig. 7. Azimuthally averaged concentric dielectric tubes as in [15] in theregiona « F < b 
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(a) (b) 
Fig. 8. Cross sectional views of anisotropically conducting tape helix supported by 
azimuthally averaged dielectric rod 
Dielectric loaded Anisotropically conducting Tape Helix model for the 
slow-wave circuit. 


Axially symmetric mode of operation 
Axially confined electron beam partially filling the tube. 
Negligible effect of transverse electric field on the electron motion 


Zero transverse speed of electrons when electron enters the interaction 
region. Richards Joe S/IEEE INAE/2016-12-13 13 
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Fig. 9. Cross sectional view of the interaction region 
Nonrelativistic operation justifying the dropping of r.f. magnetic force 
terms from the electron ballistic equation. 


The axial speed v, and the charge density o, of the entering electron 
stream remain constant with respect to the transverse co-ordinates and the 
time. 

The electron entrance speed V9 = Vy, the coldwave phase speed at the 
input signal frequency 2: 


Thermal effects are not included in this analysis. 
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B 





Axial Electric field and other field components are obtained as a 
function of electron arrival time. 


» Electron arrival time (tatz — d) vs electron entrance time 
(to at z = 0) 


> Exit electron speed (v at z = d) vs electron entrance time (tg at z — 
0) 


Parameter @ is obtained from coldwave power flow and the dc power 
required to generate the desired electron beam velocity. 


Ao, the amplitude of the r.f. input signal is determined based on a. 
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Variables rendered dimensionless 
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Dielectric 
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(a) (b) 
Fig. 8. Cross sectional views of anisotropically conducting tape helix supported 
by azimuthally averaged dielectric rod 


Electron beam radius, @ = ee Dimensional Variables Dimensionless Variables 
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a | | 

Tape helix radius, 1 — 1 radial co-ordinate 
_ time 
| b E, electron entrance ti 
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Vo 
Axial electron speed, vo(z, r, t) — 1/t;(z,r, to) (1) 
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Dimensional Variables Dimensionless Variables 
C(Z, T, tg) : time of arrival at the location t(z,r,to) = Wot(Z,T, to) 
(Z, r) of an electron with entrance time tg 


V(Z,r,tg) : axial speed at the location PUT 4) —v(ZT it, 
(Z, r) of an electron with entrance time t, 


D (Z, Tr, ts) : electron charge density Mz —p ony iç) iyii 


I(Z,Y, to) : convection current density i(z, r,t) = Vo Zo UZ,T, fo) /WoAo 


E,.(Z,7,t): for k = 1,2,3, axial, azimuthal | €,(z,7r,t) = €;,(Z,7,0)/Ay fork = 1,2,3 
and radial component of electric field vector 


H ,.(Z,7,t): for k = 1,2,3, axial, azimuthal | H,,(z,7r,t) = Z)#H;,(Z,7,0)/A, for k = 1,2,3 
and radial component of electric field vector 





Zo -free space intrinsic impedance 


Ag - axial electric field amplitude on the 


tape helix at the entrance plane 
Richards Joe S/IEEE INAE/2016-12-13 17 


Fourier series expansion of current 
density and charge density 


* qo(— v&poZo/eoAg) corresponds to both the dimensionless 
convection current density and the charge density at the entrance 
plane, z — 0. 

e tọ corresponds to the entrance time of the lf” electron which 
contributes to the time of arrival at (z,7) location at time t. 

* The current density and the charge density terms are then expanded 
in Fourier series 


i(zr,t) =Yidolte,(z.7,torlz7,t))| (2a) 
p(z.r,D 5 ; t, (z, PEZ t)) (2b) 


Jta (zr, tolz, r, t))| 
oO 


İ(Z,T,t) — qolio,a (r) + 2. (L,(Z,r) exp jmt 4 c.c.) (2c) 


m=1 
pz,7,0) = doloar) + ) (Pmlz r) expjmt + c.c.) (2d) 


m=1 
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Fourier coefficients are expanded in terms of the contribution of every 
electron’s electron arrival time 


L(Z,T) 

- (=) [Niu r, to (Z,7, 2d exp(—jmt)dt (3a) 
= (qo/2T) | a r, t9) dt, (3b) 
Py T) o 


t, (Z, FLU, t)) 


— (94/21) | t, r, ty)expC-jmt(z, r, to) dto (4b) 


T 


— = (qo/27) | > eee Bes Pi, t)) Jaime) dt (4a) 
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Approach 
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Maxwell's equations and Electron 
Ballistic equation - simplified 


A,K(T)E,, — az (Ha, + LA = —İ(Z,r,t) (5a) 


ak(r)E2, — H3, + a530,, — 0 (5b) 
a,K(r)E,, + Hz, = 0 (5c) 
aı Hı, + az (E2, + €,/r) = 0 (5d) 
4H, + €3,—a,€,, = 0 (5e) 
aı H3, — Eş —Ü or) 
C1, + (2 (Ez, T E3/T) — p(Z, r, t)/a4 (6a) 
Hı, + A2(H3.+ H3/r) = 0 (6b) 


3 i 
LL ID pe e(t,(z, ie ty)) ea EN LU ba) (7) 
with e being the electron charge and m, 1s the rest mass of the electron. 
Ay, — V9/C, 0; — (79/090), € — Ape/MeWoVo 


© {h ere | 
COE PIEDE: 
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Boundary conditions 


In the entrance plane of the interaction region, z = 0, 


t(0,r, t4) = ty (8a) 
t,(0,T,t)) — 1 (8b) 


The monochromatic signaling conditions at the tape helix 1n the entrance 
plane, with the phase factor, A( & exp(j0)), in the entrance plane is 


FE 
€,(0,1,t) = (=) exp jt 4 c.c. (8c) 


A 
€, (0,1,t) — —jf 9 exp(Jt) 4 c.c. (8d) 


Along the boundary of the electron beam radius, r — a, the continuity of the 
electric and magnetic field components are 


€,(Z,a—,t) - €y(Z7, a, t) - 0 (9a) 
Jt,(z,a—,t) - 34,(z,at,t) 2-0, k-—123 (9b) 
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At the tape helix radius, the axial and azimuthal components of the electric 
fields are continuous in (10a-b). Equations (10c-d) states the discontinuity in the 
magnetic field amounting to the surface current density of the parallel 
component. In (10e), the null electric field exists only in the tape helix region. 


E,(z,1-,t) — €,(z,14+,t) =0 (10a) 
E,(z,1-,t) — €,(z,14+,t) =0 (10b) 
(30, (z, 1—, t) — H, (z, 1+, t) ) sin Y 

t (30, z, 1—,t) — 945, (z, 1, t)) cosy — 0 (10c) 


(30; (z, 1—, t) — 96, (z, 19, t)) cos p 
— (30, (z, 1—,t) — 96, (z, 1, t)) siny — i(z r, t) (10d) 
(€4(Z, 1,t) siny — £(z, 1,t) cosylg($,z) = 0 (10e) 


where g(4, Z) is the indicator function from [13, 15, 17], which limits the 
current to the tape helix surface only. 
The boundary condition at the outer conductor,r — b —, 1s 


E,(z,b—-,t) =0, fork =1,2 
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Approach (re-view) 
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Partial Green's Function seguence from 
Fourier coefficients of Field components 


A Green’s function sequence (Gm) will be developed for the slow wave 
circuit form > 1, with the mth harmonic of the input signal frequency,mWo, 
corresponding to the mth Green’s function. 


The field components are expressed in Fourier series, 


E,(z,7,t) = [Em (Zr) expUmt)+c.c] fork =1,2,3 


M: 


E 


[Him (Z, r) exp(jmt) tc. c], (12) 


M: 


AZT) = 


m-1 
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Boundary conditions: Fourier coefficients 


By substituting the Fourier expansions of the fields (12), current density and 
charge density in the signaling conditions and boundary conditions, 


E, (0, 1) — (5 ) O1m 


Eim „(0, 1) = —] Bı (5 | Óim (1 3) 
Ej (z, a —) — Exm(Z,a +) = 0 
Hiem(Z,a —) — Hem (z,a +) = 0, k = 1,2,3 (14) 


e LE Al) — 5 
Ezm(Z, .— Ezm(Z, 14) — 
(Hım (2,1—)—H,.,.(2,1 +)) sin YU 
+(Hom(z, 1 —) — Fy 1 +)) cos V =p 
(H ne) — I +)) cos 
—( Həm (z, 1 - n Hm (Z, 1 +)) sin Y = = Lm 4 1) 


(mm pA 1) sin y — pam 2, 1) COS W)g(¢, z) = = 0 (15) 
E,,(z,b -) =0, k =1,2 (16) 
| , ifm=1 

O1m = lo. otherwise 
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E10 = H49 — E29 — Hao = 0 
H20 = 41 E30 = qo| rloa] r) + kap lr)/r|/2a2 7) 


where 
A EST 
İle) ^ otherwise 
As in [8, Il], the substitution of Fourier coefficients of the fields into the 
Maxwell’s equations (5-6) are also carried out. 
Then the non-homogeneous term i,,(Z,r) is represented as a function of the 


axial components bounded by the interval [0,d] 


Lg e > İmn(T)exp(—inkaz) for O<z<d (18) 
n=—cCo 
Ln(—n) (r) = ign (Qr) 


Lo) (7) | Tun r)cosnkgz dz (19) 
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mih harmonic of Axial Electric Field — 
Homogeneous and Particular solutions 


E, (zr) » EU? (zr) - EU) (z r) (20) 


The homogeneous component of the axial electric field 1s given as 
Ej G,r) 2 Ej. (7) exp(—7BmZ) + Ej, (r) expGf,z) (21) 


A 
f | : 
po» (z,r) = Ölm 5 -exp( ml 


ma, —— Wy) Cos Pnz if ies Vin (x, y)y dy dx (22) 


Um (x, y) = — 2 Ga (1, y) cos nkqx (23a) 
nez 
W, (r)— air e eniyi Üsr<l (23b) 
Amoo(r) / Amoo , I<r<b 
(a g^ ICA TOES n) d ; I Ban} 
a cm a) Cy (Pon a) Waa — Co (Pmn) Co (Pmn y ) (2 4) 
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Gant. y İ—— 


28 


Particular solution of the axial electric field 


E — (zr) = > Eymn(r) exp(—jnkaz) (25) 


nez 


dra 
Eimn (r) = und | | | Gant? y ) Um De y )y dy 
A -1 AL 4 
» 
— | Rinn® y ) im(x, y )y dy cos nkgx dx (2 6) 





nnn (r, y ) 


a Dm) = Ap 2 Ua 

| AR ısPmnıX D1000 I r + WC 
Tand (p r) I [0,a] (r) E (Din ) | 7 Cy (Pmn) A — Aç 1000) | 101 (r) + or) 
T in pu y Có (Pmny) aC (Pn @) Cy (Pmn®) Waa 


Amno 0 a Co (Pgh) Cy (Pmn® Waa 


Co(Pmny)l,0 sr<1 
(27) 
Co (PmnY ), 1 Lr <= b 


, G@mny) 
R mn (r, y ) — Tmn Co (nT) I [0,a] (r) YE.) (2 8) 
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Terms in the axial electrical fields 


i= Do (Dn) - go, z) : n (Pmn)Dı(Pmn) (1 v ^p1000) B 
a Co (Pmn) + AC, (Pinn) D 0 (Pmn) Cy (Pmn) ( I ^c1000) 


n us cot H Amni1 


Ania (29b) 
G5 a 
nı, D 0 (Pmn a) u D 1 [Po a) (29 C) 
we Co ( Pmn a) Cı ews a) 
Monten mn mn/Omn (30 a) 


Pmn€1Pmn)4mnoo 

Ap1000 = Cmn oi mn mno ` a (30b) 
PmnP1(Pmn)Amnoo 

Amnoo = Co(Pmn)Do(Pmnb) — Comnb)Do(Pmn) (31a) 

Amnor = Co(Pmn)D1(Pmnb) — Ci(pmab)Do(pma) | (31D) 

Amnio = C1(Pmn)Do(PmnP) — Co(Pmnb)Di(Pmn) | (310) 

Amnit = €1(Pmn)P1@mn?) ~ Gaah) Dina) (31d) 


Amnoo (r) = Co (Pnn) Do nias ) — Co (Pmnb)Do pon (32) 
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1 (29a) 


30 


Amoo = Co (Pm) Do (pmb) — Co (pmb) Do (Pm) (33) 
Amoo(r) — Co (pr) Do (pmb) — Co (pmb) Do (Dr) (34) 


Dg 5 a) PoPa) p 1] (r) (35) 


Wo(r) = 7 ee ay "nt Co(pmnr) 


where, C, and D, are common symbols for the Bessel's functions and 
modified Bessel's functions for first kind and second kind 


io, Ge 
C | Frys | v Tn J mr ‘26a 
(Pn ) J v [Dos h P < 0 l ) 
Ky p Pma?) Tinn > 0 
Dee 36D) 
Pmr) = OHE Ypma), Bins o C 
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Transverse mode number, Tmn 


ad 2 | 
Hou — co and ns — (37a) 
| 14 d. mus | 
pa = zlandpi - |t | (37b) 
2 | 
|o (à, 0€r«1 
Zar) = [n2k3 — m2azx(r)]/a3 = | mm OST ST gy 
To bern 
(2 | 
E | | te, UF 1 
t2(r) = (62, — «(r)m2a2|/a3 = | : (39) 
es. 9B 


where 68,,(> 0) is the phase shift constants for the harmonic radian 
frequencies MW, derived from the cold wave analysis of dielectric 


supported anisotropically conducting tape helix structure [13]. 
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The complete solution: Axial Electric field 
in terms of electron arrival time 


Pye yf = ala) (G JexoGict — Bız) ) + C.C. ) 


+ = -) > (>) exp(jmt) | dx 1 L Gm (Z,7T:xX,y)y dy + | Ml Zr, yay) 


f i exp(—jmt(x, y, t) )dt +c.c (40) 





where the partial Green’s function sequences G,, and R,, are given as 


Gm (z QU)X,y ) = > = (r, y ) CAP (—j nkaz ) 


nez | | E 
— Wm (r)Gmn (1, Y ) COS (BmZ) | COS (NK gx ) C29 
Ry G,rix,y) 9 9. Ros (r,y) expC-inkaz) cos(nkax) (42) 
ned 
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Electron Arrival time and Transit time 


On re-arranging and integrating the electron ballistic expression (7), the 
electron arrival time t(Z,r,tg) at any (Z,r) co-ordinate in the interaction 
region 1s obtained as (43) 


1 
rZ | (Xx 12 
toa bibendo | dx / f — 2e | €4 (s,r, t(s,r, us] 
0 | 0 ) 
forO<z<d andO<r<a_ (43) 
The electron transit time 0(Z,r,tg) is the time taken for an electron to 


reach (z,r) when it has originated the entrance plane with an entrance 
time, Lo. 


O(z,r,to) — t(z r,to) — to (44) 
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Cold wave power and input DC power 


The electron beam current of J) amps and the anode voltage Vg are computed 
accelerate the electrons from zero initial velocity at cathode to the required 
electron beam velocity vg at the interaction region's entrance plane (Z,r). This 
amounts to the DC power, Pac. The cold wave power propagating forward at 
the signal frequency wọ 1s represented as Pin[9, 12], where 


Pin = MAŞA YoPyı (45) 


with P,, being the dimensionless power computed from the cold wave fields 
from [13, 14] and Yọ is the free-space admittance (inverse of free space 
ımpedance,Z, — 1207). 


æ = 101ogio(Pin/ Pac) (46) 


The factor a in (46) provides the amplitude of the r.f. signal Ag to be chosen 
for the intended level of signal amplification as in [9, 12]. 
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Approach (re-view) 
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(the solution for the 
electron arrival time) 
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The Fourier series expansions of the axial electric field 1n (40) 1s truncated to third 
temporal harmonic (1 < m < M = 3) and 64th spacial harmonic(|n| < N = 


64) as [8, 11]. 


Operating Frequency 
Input-signal phase factor 
Tape-helix pitch angle 
Tape-helix radius 
Outer-conductor radius 


Effective dielectric 
constant of support rods 


Electron Beam current 


Electron beam radius 


Interaction region length 


Parameter, c 
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6GHz 
l 
10° 
0.01 m 
0.0224 m 
2.25 


60mA 
0.005 m 
120 
-25dB 
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The electron beam velocity required to generate the respective anode 
potential is maintained slightly above the cold-wave velocity v,, 
obtained from the dispersion characteristics plot[14, 15]. 


With v, — 0.13579066 * c where c is the velocity of light, the electron 
beam velocity at the entrance plane 1s chosen slightly greater than v, as 
Vo = 0.1394407144 x» c. 


The normalized propagation phase constants from the dispersion plot 
of anisotropically conducting tape-helix supported by dielectric rods [13, 
14] are 54 — 1, 6, — 1.6127704 and 4, — 2.225748. 


On solving the electron arrival time (43) by successive approximation 
technique, the electron arrival time at the end of interaction length (d) 1s 
plotted. 
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Electron Transit Time 


The media below displays an example of electrons entering at 
entrance plane z — 0 with different entrance times, take 
different transit times to reach the exit plane z — d. 


t- t0 t(d.r.t0) 


| 
| 
| i 


AKI ui RI 
| 


z-0 2*pi z=d 








Fig. 9 The electrons arriving at different entrance times, travel at different velocities to 
reach the exit plane at different arrival times 
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Electron transit time, O(d,r,t,) 


118 





116 | 
—T -7/2 0 TU A T 


Electron entrance time, to 
Fig. 10 Electron transit time at the exit plane, plotted against electron entrance time for the electrons 


at the beam center(r — 0), half way point, (r — a/2) and at the edge of the electron beam(r — p 
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Electron arrival time t(d,r, tg) 


Electron arrival time, Kdır,t, ) 
oc 


1146> - 





114 | | 
—T -7/2 0 TU 2. T 


Electron entrance time, ^ 


Fig. 11. Electron arrival time at exit plane, plotted against electron entrance time for the electrons at the 


beam center(r = 0), half way point, (r = a/2) and at the edge of the electron beam(r = a) 
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Radial position, r 


Electron bunching with a = —25dB 


e Electron bunching along the interaction length is plotted below, for a = 
— 25dB. 
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Fig. 12 Numerically computed Electron arrival time along z at various 
radial locations are plotted 
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1.04 


Electron exit velocity, v(d.r.t) 


0.94 





—T —7/2 0 T/ 2 n 
Electron entrance time, t 
Fig. 13. Electron exit velocity at the exit plane, plotted against electron entrance time for the electrons 


at the beam center(r = 0), half way point, (r = a/2) and at the edge of the electron beam(r = eh 
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Lİ A rigorous large signal field analysis was performed on the 
practical case of linear beam travelling wave tube with tape 
helix slow wave structure supported by dielectric rods. 


Lİ The transit time and the exit velocity at the exit plane of the 
interaction region were computed and plotted. 


L] The numerical computation of the axial electric field, 
current gain, power gain and the conversion efficiency are 
in progress. 


Lİ In the future, the large signal field analysis will be performed 
in a 3D simulation software. 


Richards Joe S/IEEE INAE/2016-12-13 46 





The manuscript presents the computation of the electron transit time and 
the electron exit velocity in a dielectric supported helix traveling wave tube 
amplifiers (TWTA) by solving the non-homogeneous boundary value 
problem. The following suggestion can be incorporated to bring out the 
novelty of the presented work. 


1. The large signal analysis results of tape helix model analyzed in this 
paper could be compared with the results obtained from the existing 
analysis of sheath helix model. 


Since the numerical computations include rigorous field analysis, the 
analysis results of sheath helix model could not completed in the short 
duration of time. However, the suggestion will be considered for the future 
works. 


2. Maintain consistency in referencing symbols throughout the manuscript 


The notations of variables were updated and corrected in the 
presentation. 
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